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Abstract 

The recent TRIUMF experiment for (jTp — * nVf/y gave a surprising result 
that the induced pseudoscalar coupling constant gp was larger than the value 
obtained from fj,~p — > nv^ experiment as much as 44 %. Reexamining contri- 
bution of the axial vector current in electromagnetic interaction, we found an 
additional term to the matrix element which was used to extract the gp value 
from the measured photon energy spectrum. This additional term, which 
plays a key role to restore the reliability of gp{— 0.88m^) = 6.77^(0), is 
shown to affect the gp quenching problems in nucleus. 
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I. INTRODUCTION 



The matrix element of vector and axial vector currents are generally given as 

(N(p')\V»(0)\N(p)} = u(p')[G v (q 2 )Y + ^SV + G M (q 2 )a^q u ]fu(p) 
(N(p')\AZ(0)\N(p)) = u(p')[G A (q 2 )r + ^fgV + G T {q 2 )a^q v ] l5 ^ u{p) , 



where G A (0) = g A (0), G M (0) = g M (0), G v (0) = g v (0) and G P (q 2 ) = (0g P (q 2 ) with the 
nucleon and muon masses, m and m M . r a is the isospin operator. G$ and Gt belong to the 
second class current which has a different G-parity from the first class current, and they are 
assumed to be absent from the muon capture to be discussed in this paper. On the basis of 
PCAC (Partially Conserved Axial Current), the induced pseudoscalar coupling constant is 
calculated as 



g P (-0.88m 2 ) = - ™ n - g A (0) = 6.77^(0). (2) 



i 2 ' 

'" m2+0.88m2' 



This value is confirmed by an experiment of the ordinary muon capture (OMC) on a proton, 
(Tp ^nv^ 0. 

However, in order to obtain more precise data, the TRIUMF group measured recently 
the photon energy spectrum of the radiative muon capture (RMC) on a proton, fi~p —>■ nv^ 
and extracted a surprising result 

g P = g P {-0.8Sml)/g A {0) = 9.8 ± 0.7 ± 0.3 . (3) 

It exceeds the value obtained from OMC as much as 44%. This discrepancy is serious because 
the theoretical value of gp is predicted in a fundamental manner based on PCAC and agrees 
with the OMC value. As long as PCAC is assumed to be creditable, a doubt may be cast 
on the result of TRIUMF experiment. Recent calculations 0|| by chiral perturbation also 
says such a doubt. However, in order to solve this puzzle, one has to reexamine carefully 
the Beder- Fearing formula which is a phenomenological model, used to extract the gp 
value from the measured RMC spectrum. 



In finite nuclei, through the theoretical analyses of OMC experimental results, it is 
already reported H that the gp value is quenched in medium-heavy and heavy nuclei while 
it is enhanced in light nuclei. Since these analyses are carried out before the recent TRIUMF 
experiment one needs to reconsider those analyses from another viewpoint. 

In this paper we present more successful analysis for the recent TRIUMF data and show 
some progressive results for gp quenching problems in nucleus by applying our results on 
proton to nuclear matter. 



II. BASIC FORMULAE 

We start from the ordinary linear-cx model ; 

C = n^d, - g(a + it ■ 7f T5 )]* + l -[{d^f + («V) 2 ] + + ° 2 ) ~ + v 2 f (4) 

, which gives the following axial current 

K = %M75y* + ?r a <V - ad,ix a . (5) 

By the spontaneous breakdown of chiral symmetry, a field is shifted to a' = a — Co with 
°o = /tt- Consequently, the pion appears as Nambu-Goldstone boson. The PCAC can be 
satisfied by the additional inclusion of the explicit chiral symmetry breaking term as well 
known. 

But the axial current 

Al = *7 M 75y* - UdX (6) 



gives gA = 1 in the tree approximation. Following the recipe of Akhmedov [|T^] to cure this 
problem, we add chiral invariant lagrangian C\ to Co, 

Ci = C[* T ^*0? x c^tt) + * 7AtT 5^(7fa M (7 - adrf)\ , (7) 

where arbitrary parameter C is determined so that the axial current pertinent to nucleons 
in C = C + C\ 



(Ar) ^ = ^75^[l + C 2 (7r 2 + a 2 )] (8) 

should satisfy ^ A a = Qa^IixI^ 1 ^ ^ with = 1.26. The Goldberger-Treiman relation then 
is satisfied exactly. As a consequence, ^ A^ includes the contribution not only from the 
nucleon but also from the n — N interactions. 

Now, the axial vector current consists of the nucleon and pion sectors as 

A£(x) = WAZ(x) + WA£(x) (9) 

= wa£(x) + f n dy a (x) , 

where f n is the pion decay constant. <p a (x) is the pion field. 

To describe RMC, we need a radiative axial current, which is used to obtain the transition 
amplitude of RMC by coupling to the weak current of lepton line. Three different methods 
are considered in order to construct a radiative axial current. The 1st method M is to 
start from the above lagrangians C = Co + £i using covariant derivative, by which we 
introduce a photon field in U(l) gauge invariant way. The outcoming lagrangian gives 
a radiative axial current, which characteristic is its non-conservation through the explicit 
chiral symmetry breaking due to the electromagnetic interaction. The second method || 



is to use the extended Euler equation Jl(J for the lagrangian C = C + L\. The third is 
to make it directly from the above axial currents, eq.(9), by exploiting it minimal coupling 
scheme to the momenta of relevant particles. Here we follow the third one. Of course the 
final radiative axial currents from these three different methods turned out to be equivalent. 
Let us begin from the divergences of ( N ^A%(x) and ^'A^(x), 

d^A^x) = d»[g A *(x)y»y 5 ^(x)} = , (10) 

d^A^x) = f n d 2 Mx) = -MmlMx) + Ja] , (11) 
where we used pion field equation from the above lagrangians 

(9 2 + m 2 )0 a = - Jf . (12) 



The quantity denotes the pion source term. Therefore, the divergence of total axial 
currents is given in the following way 

d.A^x) = d^A^x) + d^A^x) (13) 
= -M™&a(x) + J?) + f 1[ J? 
= -f n ml(j) a (x) . 

This is PCAC. By eqs. (10) and (12), one can obtain 

*•<*> = "T^W^M- (14) 

Substitution of eq.(14) into eq.(9) yields 

KV) = {N) K{x) - w ^{idY[d^A^{x)]. (15) 

In order to clarify a role of the axial current in description of the radiation process, we adopt 
the minimal coupling prescription, i.e. d\ — > d\ — ieA\ and — > — eA^. This procedure 
leads to 



2 (idr[ex {N) A x a (x)] + ^-^[e^Al{x)\ , (16) 

where the potential A\ is replaced by the photon polarization vector e\. Notice here that 
dx^A^(x) = g A m^>(x)i'yz ) T a ^>(x). The last term in eq.(16) can be neglected because it 
appears at 0(e 2 ) order. Thus we can express the axial current in the radiative processes in 
the following realistic form 

A£{x) = ^{x)[g Al ^ + ^V 75 - ^^75]^) (17) 
" n ' [(j2) q^(x)e al "l^(x)]. 



2mm (J 2 



The fourth term, which corresponds to "Seagull term" , is missing |5| in the previous calcu- 
lations g,[7|. 



Following the Fearing's formulation and notation || for the diagrams given in ref. 
one can evaluate the relativistic amplitude of RMC on a proton as 



M fi = M a + M b + M c + M d + M e + AM e (18) 



with 



M a = -e a u n T s (Q)u p ■ u uls {l - l5 f |±^ 7% . (19) 

—2k ■ /i 

M 6 = e a L s u n T\K)tll±^hL^ _ , 



M c = e a L s u n {-in n - hp) — — r (K)u p , 



2Q a + k a g P {K 



'2\ 



M d = -e a L 5 u n (^ — '__ ^^\- J K s >y 5 )u 

Q2 _ m 2 m ^ 
nj T - f^M Sa , 9p(Q 2 ) Sa\ 

M e = e a L s u n {—a + ^ 5 g )u p , 

2m m p 



AM e = -e a L 5 u n (^^Q s ^ a )u p , 
2mm 



f 



where 

T\q) = g vl s + f^v s % + g A i^ + ^V 7 5 , (20) 

Ls = u u , ~fs{^ — Ihju^K = n — p + k and Q = n — p with momenta of neutron, proton 
and photon, n,p and k, respectively. And m ~ m p ~ m n . Other constants are taken as 
gv = 1-0, 5u = — 1.25, #m = 3.71,/%, = 1.79 and n n = —1.91 ||. M e term is originated 
from the third term in eq.(17) and AM e term comes from the fourth term. But the latter, 
AM e , is missing in the paper by Fearing ||[7|]- Accordingly, this term was not included in 
the previous procedure of extracting gp value from the experimental RMC photon energy 
spectrum 0. 

The above transition amplitude can be also understood in terms of pseudo vector (PV) 
coupling scheme between nucleons and virtual pion, through which external axial current 
interacts with nucleons. Moreover in nuclear matter, this PV coupling type is preferred 
rather than PS coupling type because the former is consistent with PCAC while the latter 
contradicts to PCAC in nuclear matter [|l^ . 



The RMC transition rate is given by 

dT RMC aG 2 \^\ 2 m N f 1 kE 2 1 „ 2 

where a is the fine structure constant, G is the standard weak coupling constant, y = 
k ■ v, k max = - m 2 n )/2W , E v = W (k max - k)/[W - k{l - v)), W = m p + m n - 
(muon binding energy) and \<p^\ 2 is the absolute square of muon wave function averaged 
over the proton which is taken as a point Coulomb. In order to compare to the experimental 
results, we take the following steps. For liquid hydrogen target, muon capture is dominated 
through the ortho and para pfip molecular states |p|p!3|. Since these molecular states can be 



attributed to the combinations of hyperfine states of fip atomic states |13[ i.e. single and 
triplet states, we decompose the statistical spin mixture \ J2 sp ins l^/«| 2 into such hyperfine 
states by reducing 4x4 matrix elements to 2 x 2 spin matrix elements. At this step, we 
confirmed that when the AM e term was not included, eq.(21) reproduced the curves given 
in ref. 17H. 



For the description of RMC in nuclear matter, we follow the Fearing's paper [T4 



adopted the relativistic mean field theory |T5] where the nucleons are treated as free Dirac 
particles with effective mass due to the scalar and vector potential. Then the nucleus are 
the Fermi gas. Our RMC capture rate in nuclear matter is given in the following way 

^ NM aG 2 3m P m N 2 f k ^ [i cose k) min r k F 



T rmc = 4 2 4 , 3 1^1 / dk . . . d ^s6 k / dp dn 

f 2n ,, E v np 2 k 1 ^ l9 . . 

/ dfa - E \MM 2 > (22) 

Jo E p E n \n + p\ A spins 

where the integration intervals of nucleon momenta, which comes from the Pauli blocking in 
nuclear matter, are calculated in detailed kinematics and kp is the Fermi momentum. For 
finite nuclei, we need to know the corresponding muon wave functions and kp values, but 
which depends on the model. We already suggested a model for this purpose, but will 
be skipped here and concentrate on the case of nuclear matter. 
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III. RESULTS AND DISCUSSIONS 



Our results for RMC on proton are shown in Fig.l. The solid curve is the spectrum 
obtained in ref. 2., i.e. the result without AM e term for gp = 9.8. On the other hand, 
the dotted curve is calculated without AM e term for gp = 6.77. This curve is obviously 
much lower than the solid curve. When AM e term is taken into account for gp = 6.77, we 
obtain the dashed curve which is very close to the solid curve for the energy spectrum on 
k > 60MeV. The minor discrepancy may be due to the neglect of higher order contributions 
and other degree of freedom such as A. Our result shows that AM e term restores the credit 
of gp = 6.77. 

The number of RMC photons observed for k > 60MeV is 279 ± 26 and the number 
of those from the solid curve is 299, while our result obtained by integrating the dotted 
curve spectrum is 273. Since the contribution of A degree of freedom is known to be a few 
percent 0, it is not included in the present calculation. Vector mesons such as p and u are 
also turned out to have very small contributions in this calculation. Higher order terms are 
pointed out to be insignificant ||. 

It is confirmed that the matrix elements upto the M e term in eq.(19) satisfy gauge 
invariance. However, it is broken when the AM e term and A degree of freedom are included. 
As far as diagram method is adopted, the gauge invariance is more or less broken because a 
series of diagrams will be cut somewhere. In order to estimate the rate of gauge invariance 
broken by the AM e term, we evaluated the spectrum with AM e alone. The result is shown 
by a dot-dashed curve in Fig.l. The gauge invariance breaking is not so large as expected. 
However, in the limit m M — > 0, the gauge invariance is restored. 

In spite of such a burden of gauge invariance, our present calculation shows that gp = 
6.77 is reasonable for both OMC and RMC on a proton. 

Figure 2 shows the photon energy spectrum in nuclear matter, which is presented as the 
ratio of RMC to OMC in order to reduce the uncertainty from the nuclear structure. We 
compared our amplitude to Fearing's analysis, which is PS coupling and has been used to 



8 



extract the gp from the experimental data. At the same gp, our ratio R is higher than the 
PS coupling scheme. It means that gp value to fit some experimental data becomes lower in 
our PV scheme. Therefore gp quenching rate is smaller than the usual PS coupling scheme. 
These behaviour are nearly independent of effective nucleon mass in nuclear matter (see the 
solid, dotted and dashed lines). 

Figure 3 shows another interesting results for gp quenching in finite nuclei. The larger 
which may mean the heavier nuclei, the smaller becomes the ratio R. As a result, gp 
quenching may be larger in heavier nuclei. But in relatively lower kp region just reversed 
results are shown. Consequently gp may be enhanced in lower kp region. Since we have 
to integrate the nucleon's Fermi motion from the possible lowest momentum up to the 
Fermi momentum, the mechanism in lower kp region could play a role of compensating gp 
quenching in larger kp region. Therefore this could be an indication for the gp enhancement 
in light nuclei. 
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Figure Captions 

Figure 1. Photon energy spectrum for triplet states in RMC on proton. The solid curve 
is deduced without AM e term for gp = 9.8, whose result corresponds to the experimental 
results in ref.2. The dotted curve is obtained without AM e term for g P = 6.77. The dashed 
curve is with AM e for g P = 6.77. The dot-dashed curve is calculated with AM e term alone 
for gp = 6.77. 

Figure 2. The photon energy spectrum for the ratio of RMC and OMC in nuclear matter. 
The thick curves are results from our transition amplitudes, but thin curves are Fearing's 
amplitudes. The solid curves are for M* = M, the dotted curves are for M* = 0.57M and 
the dashed curves are for M* = 0.7 M. 

Figure 3. The ratio of RMC and OMC versus fermi momentum hp. The dot-dashed 
curve (fc 7 = 60 MeV) and dotted curve (/c 7 = 80 MeV) come from our amplitude, while the 
solid and long dashed curves are from Fearing's, respectively from k 7 = 60 and 80 MeV. 
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